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Constraining the population of cosmic ray protons in cooling 
flow clusters with y-ray and radio observations: 
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Abstract. We wish to constrain the cosmic-ray proton (CRp) population in galaxy clusters. By hadronic interactions with 
the thermal gas of the intra-cluster medium (ICM), the CRp produce y-rays for which we develop an analytic formalism to 
deduce their spectral distribution. Assuming the CRp-to-thermal energy density ratio Xqr p and the CRp spectral index to be 
spatially constant, we derive an analytic relation between the y-ray and bolometric X-ray fluxes, T y and F x . Based on our 
relation, we compile a sample of suitable clusters which are promising candidates for future detection of y-rays resulting from 
hadronic CRp interactions. Comparing to EGRET upper limits, we constrain the CRp population in the cooling flow clusters 
Perseus and Virgo to X CRv < 20%. Assuming a plausible value for the CRp diffusion coefficient k, we find the central CRp 
injection luminosity of M 87 to be limited to 10 43 erg s -1 a-/(10 29 cm 2 s~'). The synchrotron emission from secondary electrons 
generated in CRp hadronic interactions allows even tighter limits to be placed on the CRp population using radio observations. 
We obtain excellent agreement between the observed and theoretical radio brightness profiles for Perseus, but not for Coma 
without a radially increasing CRp-to-thermal energy density profile. Since the CRp and magnetic energy densities necessary 
to reproduce the observed radio flux are very plausible, we propose synchrotron emission from secondary electrons as an 
attractive explanation of the radio mini-halos found in cooling flow clusters. This model can be tested with future sensitive 
y-ray observations of the accompanying 7r°-decays. We identify Perseus (A 426), Virgo, Ophiuchus, and Coma (A 1656) as the 
most promising candidate clusters for such observations. 

Key words, galaxies: cooling flows - galaxies: cluster: general - galaxies: cluster: individual: Perseus (A426) - intergalactic 
medium - cosmic rays - radiation mechanisms: non-thermal 



£> . 1. Introduction lBrunettil2002l) which produce radio signatures that differ mor- 

phologically as well as spectrally: 

Cooling flows are regions where the influence of non-thermal 



C3 ■intra-cluster medium (ICM) components such as magnetic 1- Shock acceleration: Natural acceleration mechanisms pro- 
fields and cosmic rays may be strongest within a galaxy cluster viding relativistic partic les are strong structure formatio n 
owing to strong observed magnetic fields, central active galax- and merger shocks (e.g jHarris et alJll980tlSarazinlll999l) . 
ies, and increasing non-thermal-to-thermal energy ratio due to Detailed studies h ave been undertaken on shocks of cos- 
rapid thermal cooling processes. They are also regions where mological scales JMiniati et alJ l200fj l^kizawa & Naitd 
such components are best detectable due to the high gas den- 2000[). Fermi I acceleration processes of CRe at these 
sity which allows for secondary particle production in hadronic shock fronts produce la rge scale extended p eripheral ra- 
interactions of cosmic ray nuclei with the ambient gas. By the dio relics as proposed by En Blin et al]dl998al). For instance 
term cooling flow we do not rely on specific models but only two prominent relics in Abell 3667 jRottgering et al J 19971) 
on observed properties such as declining temperature gradients were s uccessfully reproduced in a simulation of this pro- 



and enhanced electron density profiles towards the center of the cess bv lRoettiger etal.l ( fl99 

cluster. 2. Reaccelerated electrons: Secondly, reacceleration pro- 

XT , , .. . ,. t . , , .. , cesses of mildly relativistic CRe (y ^ 100 - 300) be- 

Non-thermal relativistic particle populations such as cos- J v/ ' 

i » /rin \ a . /rir> \ u ■ ■ . j ina injected over cosmological timescales into the ICM by 

mic ray electrons (CRe) and protons (CRp) can be injected b J b 1 

. T< ^ A . . , , ,-rc t , t „ . sources like radio galaxies, supernova remnants, merger 

into the ICM mainly by three different processes (following b ' * ' b 

shocks, and galactic winds can provide an efficient sup- 

ply of highly-energetic CRe. Owing to their long lifetimes 

Send offprint requests to: Christoph Pfrommer, of a few times 10 9 years these mi ldly relativisti c CRe can 

e-mail: p£rommer@mpa-garching . mpg.de accumulate within the ICM (see Sarazin 2002, and ref- 
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erences therein), until they experience continuous in-situ 
acceleration either via shock acceleration or resonant pitch 
angle scattering by turbulent Alven wave s as originally pro- 
posed by Ijaffd d 19771), reconsidered bv ISchlickeiser et alJ 
dl987l) . and lately by bhno et all d2002l) . These accelera- 
tion processes of CRe possi bly yield exten ded radio halos 
centered on the cluster ( Brun etti et alJ20 01 ) while there are 
also suggestions that radio mini-halos wi thin a cooling flow 
cluster originate from these processes ddittietalJ 12001 . 
There is also evidence that reacceleration processes act- 
ing on fossil radio plasma produces small filamentary radio 
relics at the cluster periphery, so-called revived ra dio ghosts 
fenfilin & Gop al-Krishnal l200lt lEnfilin & Briiggenl 12002) 
presumably by adiabatic compression in shock waves. 
3. Particles of hadronic origin: Eventually, CRp can in- 
teract hadronically with the thermal ambient gas produc- 
ing secondary electrons, neutrinos, and y-rays in inelas- 
tic collisions taking place throughout the cluster volume 
which would generate radio halos through synchrotron 
emiss ion (first pointed out by iDennisonl Il980t IVe strand 
1982). In the ICM the CRp have lifetimes of the or- 
der of the Hubble time JVolk et al.l Il996l lEnfilin et alJ 
ll997tlB~erezinskv et alll997l) 7 long enough to diffuse away 
from the production site and to maintain a distribution 
over the cluster volum e. This process w a s rec onsidered 
in more detail by iBlasi & Colafrancescol d 19991) and by 
iDolag & Enfilinl J2000l) . the latter authors using numeri- 
cal hydro-dynamical simulations including magnetic fields. 
Recently, Mini ati et al 1 ll2001bl) have performed cosmologi- 
cal simulations of cluster formation including injection pro- 
cesses of primary CRp. These authors conclude that under 
certain conditions extended diffuse radio emission could be 
due to hadronically produced CRe. However, there are also 
claims that extended radio halos cannot be generated by 
secondary electrons due to the morphological steepness of 
predicted radio brightness profiles in contrast to observa- 
tions (Brunetti 2002). Besides constraining the population 
of CRp in the ICM, this work will present arguments for the 
hadronic origin of radio mini-halos or a substantial contri- 
bution of secondary electrons to these mini-halos. We fur- 
ther perform a parameter study which shows that the large 
cluster radio halos could be also of hadronic origin, pro- 
vided the CRp-to-thermal energy density profile is radially 
increasing. 

It is very difficult to distinguish between contributions of 
these three populations of cosmic ray (CR) particles to non- 
thermal particle populations, especially if all of them account 
for injection of cosmic rays into the ICM in different strength 
depending on underlying governing physical processes and pa- 
rameters. The hadronically produced CRe may be reaccelerated 
by shocks or cluster turbulence and therefore mix up the differ- 
ent CRe populations. 

Radio observations of the radio halo in the Coma cluster 
find a strong steepenin g of the synchr otron spectrum with in- 
creasing radius dGiovannini et aljfl993l) . This b ehavior is ex- 
pected for a reaccelerated population of CRe dBrunetti et alJ 
1999J l200ll) . There is also a report of radial spectral steep- 



e ning in the case of the radio mini-halo of Perseus according 
to Siibring lll993l) . This, however, could easily be an observa- 
tional artifact owing to a poor signal-to-noise ratio in the outer 
core parts of the cluster in combination with the ambiguity of 
determining the large scale Fourier components owing to the 
nonuniform coverage of the Fourier plane and missing short- 
baseline information: the so-called "missing zero spacing"- 
problem of interferometric radio observations. By comparing 
the spectral index distribution of the three radio maps (92 cm, 
49 cm, and 21 cm), there seems to be likewise a possibility of 
radial spectral flattening depending on the chosen radial direc- 
tion. The hadronic electron model does not necessarily produce 
the radial spectral steepening without fine-tuning. 

Assumptions: The purpose of this work is to provide con- 
ceptually simple analytic instruments for describing the spec- 
tral signatures in radio, X-rays, and y-rays resulting from in- 
elastic cosmic ray ion collisions. It is especially important to 
constrain the population of CRp within clusters of galaxies in 
order to understand the governing physical processes of these 
objects and the important theoretical implications for the non- 
thermal content of the ICM, i.e. if non-thermal CR pressure 
plays an important role in supporting the intra-cluster ionized 
gas dEnfilinet al.ll997l) . The assumptions of our models are: 

- CRe are taken to originate from hadronic interactions of 
CRp with thermal ambient protons of the ICM and the CRp 
population is described by a power-law distribution in mo- 
mentum. The origin of this population is not specified here, 
but CRp may be accelerate d by shock waves of cluste r 
mergers, accretion shocks ( Cola francesco & Blasil fl 998). 
or injected from radio galaxies into the ICM ( Valtaoia 1984; 
EnBli n et al.lll997t IBlasi & Colafrancescoll 19991) or res ult 
from supernova driven galactic winds dVolk et al.ll99rj|) . 

- In our isobaric model, the energy density of CRp is as- 
sumed to be proportional to the thermal energy density of 
the ICM. In our scenario of adiabatic compression of CRp 
during the formation of the cooling flow this proportion- 
ality is imposed prior to the transition. This assumption is 
reasonable if the thermal electron population and the CRp 
were energized by the same shock wave assuming that there 
is a constant fraction of energy going into the CRp popula- 
tion by such an acceleration process. As a third model we 
take a single central point source injecting the CRp which 
result s in a very peaked CRp profile (compare lBlasilll999t 
IBlasi & Colafr ancescoll999l) . 

- The CRp spectral index is assumed to be independent of 
position and therefore constant over the cluster volume. 
In some sense this represents an oversimplification which 
could be abandoned in order to reproduce some specific ob- 
servational results, which however would be questionable 
without understanding the underlying physical processes. 

- The electron density and temperature profiles of the ICM 
are assumed to be spherically symmetric and were taken 
from the literature. This assumption is justified in the case 
of y-rays resulting from neutral pion decay because we use 
cluster volume averaged spectra in order to compare to ob- 
servation, and is not severe in the case of radio emission, 
since the profiles are obtained from deprojected X-ray data. 
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The magnetic field configuration is assumed to be spheri- 
cally symmetric on cluster core scales and follows the elec- 
tron density with a power- law index as a free par ameter 
within the suggested range (Dolag et a ll200ll lT999). 

- No reacceleration or diffusion process of CRe is taken 
into account in calculating the synchrotron and the inverse 
Compton (IC) emission. Therefore we provide conservative 
estimates for the flux. 

- The radio spectrum is taken to be quasi-stationary owing 
to the short electron cooling time which establishes a sta- 
tionary CRe population on very short timescales. There is 
a one-to-one correspondence between the CRp power law 
index and that of the CRe population which is in addition 
determined by radiative synchrotron losses and IC cooling. 

The paper is organized in two main parts: It starts with 
theoretic modeling of y-ray spectra resulting from hadronic 
CRp-p interactions and presents models of synchrotron and 
IC radiation emitted by secondary electrons. The second part 
discusses the astrophysical application of this formalism to a 
nearby cluster sample including cooling flow clusters. After 
modeling the spatial distribution of CRp within cooling flow 
clusters we constrain this population by comparing to EGRET 
upper limits. We furthermore obtain limits on the CRp pop- 
ulation by the morphology of radio brightness profiles in the 
case of Perseus and Coma. Throughout this paper we assume 
the standard ACDM cosmology with Qm = 0.3, £2a = 0.7, 
and Ho - lOhio km s~ Mpc~', where h-ia indicates the scaling 
with Ho- 



2. Theoretic modeling of multi-frequency 
signatures resulting from hadronic CRp 
interactions 

In order to study non-thermal emission from clusters we model 
the IC and synchrotron radiation of secondary CRe produced 
in inelastic collisions by CRp scattering off thermal nuclei as 
well as the y-ray spectrum produced by decaying pions be- 
ing produced by these CRp-p collisions. After introducing our 
definitions (Sect. 12. 1> . we develop an analytic formalism de- 
scribing the decay of secondary neutral pions into two high- 
energy y-rays (Sect. l2~2l . Section l2.2.2l uses the analytical fire- 
ball model for inelastic CRp interactions with nuclei of the in- 
tergalactic mediu m in the high-energy regime of C Rp (E p » 
m p c 2 ), following Mannheim & SchlickeiseJ 111994) . Based on 
that we develop in Sect. 12.2.31 an analytic formula describing 
the y-ray spectrum by parameterizing important effects near 
th e pion thr eshold using an approximate descrip tion developed 
bv lDermenJl986alb J ). wh ich combines is obaric t Stecker 1970) 
and scaling models (Badhwar et al. 1977; Stephens & Badhwar 
119811) of the hadronic reaction. Using this formalism, an an- 
alytic Ty-Fx scaling relation is derived in the framework of 



a simple scenario of spatial distribution of CRp (Sect. 12.31 . 
Finally, Sect. I2.4l deals with radio and X-ray emission of sec- 
ondary electrons being produced by decaying charged pions. 



2.1. Definitions 

Throughout the paper we use the following definitions for the 
differential source function q{r, E), the emissivity j(r, E) and 
the volume integrated quantities, respectively: 



q(r,E) = 
Q(E) = 



AN 



dfdVdE ' 
J dVq(r,E), 



j(r,E) = Eq(r,E), (1) 
J(E) = E Q(E) , (2) 



where N denotes the integrated number of particles. From 
the source function the integrated number density production 
rate of particles A.(r), the number of particles produced per 
unit time interval within a certain volume, and the particle 
flux T can be derived. The definitions of the energy weighted 
quantities are denoted on the right hand side, respectively, 



A(r) = J dEq(r,E), 
L = J dVA(r), 

AnD 2 



Mr) 
L 

F = 



= J dEEq(r,E) : 

= J dVA(r), 
L 



4ttD 2 



(3) 
(4) 
(5) 



2.2. y-ray spectrum from hadronic CRp interactions 

2.2.1 . Cosmic ray proton population 

The differential number density distribution of a CRp popula- 
tion can be described by a power-law in momentum p p , 

f p (r,p p )dp v dV = ncRp(r) (^f (^) dV , (6) 

where the tilde indicates that «cRp is not a real CRp number 
density while it exhibits those dimensions. We choose the nor- 
malization ncRp(f) in such a way that the kinetic CRp energy 
density BcR P (r) is proportional to the thermal energy density 
£thWof thelCM, 

ecRpW = X C R P (r) e±(r) = I dp f p (r, p p ) E m (p p ) (7) 

Jo 

= h CRp (r)m p c 2 t m p c 2 ^ ap ( a p - 2 3-or p \ 
2(a p - 1) \ GeV / \ 2 ' 2 / ' 

The kinetic energy of CRp E^ n and the thermal energy density 
of the ICM Sfh are given by 



Ekin(pp) = ^p 2 p c 2 + m 2 p 



c — m D c 



s±(r) = -d e n e (r)kT e (r), 

u j 1 1 ~ * Xlie 

where a e = 1 H ■ 

1 _ 5^ He 



(9) 
(10) 

(11) 



counts the number of particles per electron in the ICM us- 
ing the primordial 4 He mass fraction Xn e = 0.24, and &(a, b) 
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Fig. 1. The ratio of CRp energy densities scRp(Pmin, <* p ) with 
and without a lower cutoff p m [ n in the CRp number density dis- 
tribution function as a function of p m \ n for different values of 
the CRp spectral index a v (see eq. dl5l l). For CRp the kinemati- 
cally allowed threshold in order to produce 7r°-mesons hadron- 
ically is given by p^ = 0.78 GeV c _1 . 



denotes the beta-function lAbramowitz & Stegunl[l9 65). The 
functional dependence of the CRp scaling parameter XcR P (r) is 
a priori unknown. In order to draw astrophysical conclusions 
for the CRp population in clusters of galaxies, we adopt three 
different models for the spatial distribution of CRp later on in 
Sect.rO 

In contrast to relativistic electrons which loose their energy 
on relatively short time scales compared to the Hubble time 
through synchrotron emission in cluster magnetic fields and 
IC scattering with photons of the microwave background, the 
dominant energy loss mechanisms of CR p are electronic ex- 
citati ons in the pla sma (EnB lin et aljfl9 97). defining a cooling 
time (Gould 1972) 



df 



1 

r P 

m e c 3 m p J3 p y p 



4 71 e 4 « e 



In 


'2y p m e c 2 /3p 


#1 






2 



(12) 



(13) 



where fi p c denotes the velocity of the proton, y p its relativistic 
Lorentz factor, and cj p \ = (47re 2 n e /m e ) 1 ^ 2 the plasma frequency. 
Inserting typical values for cooling flows yields a lower cutoff 
on the CRp momentum 



Gyr/\10 2 cm 



GeVc" 



In general, this gives rise to a spatially dependent cutoff of the 
CRp momentum which increases with time. 

In order not to rely on too many assumptions, we do not 
impose a specific momentum cutoff which is possible since 
the spectral index a p varies in our model in between 2 and 3. 
Instead, we quantify the influence of a lower cutoff /? m ; n on the 
population of CRp by taking the ratio of CRp energy densi- 
ties ecRp(/?min, ffp) with and without a lower cutoff. This ratio 



as shown in Fig. [2 can be written using the definition for the 
normalized lower CRp momentum cutoff p = 



eC Rp(P,«p) s,(^,^)+2/5 i - P (vrT^-i) 



scr p (0, a p ) 



(15) 



where B r (a,b) denotes the incomplete beta-function 
jAbramowitz & Stegunl 19651 with x = (1 + p 2 Y l . Combining 
Fig. [0 and eq. J 141 demonstrates the small influence of 
Coulomb cooling to the CRp energy density within cooling 
flows. 

2.2.2. Fireball model 

The CRp interact hadronically with the thermal background 
gas and produce pions with relative multiplicities ^ = ^„± 
according to isospin symmetry and assumin g thermal eq uilib- 
rium of the pion cloud in the center of mass jFerrnill9 50). The 
charged pions decay into secondary electrons (and neutrinos) 
and the neutral pions into y-rays: 



2y. 



+ V e /V e + V M + Vfj 



Only CRp above the kinematic threshold p t hi = 0.78 GeV c~' 
are able to produce pions hadronically and are therefore visible 
through their decay products in both the y-ray and radio bands 
via radiative processes. Only the CRp population above this 
threshold is constrained by this work while the lower energy 
part of this population in general can not be limited by only 
considering hadronic interactions. 

In the high-energy limit for CRp (E p » m p c 2 ) the pion 
source function resultin g from hadronic CRp-p interaction s can 
be calculated following Man nheim & SchlickeiseJ (fl994) to be 



q n o{r,E n o)dE n odV ~ 2 3 cr pp c n N (r) 
<6E„o 



«CR P (r) 



GeV 



GeV 

d£j> dV, 



(16) 



where a y = 4/3 (a ? - 1/2), <x pp = 32 mbarn is the inelastic p-p 
cross section, and n N (r) = d Vdr n e (r) = « e (r)/(l - jX He ) is the 
target nucleon density in the ICM. The 7r°-decay induced om- 
nidirectional (i.e. integrated over 4 n solid angle) differential y- 
ray source function can be calculated in this energy regime as- 
suming the decay products are distributed isotropically in their 
rest frame, yielding 



(14) q y (r,Ey) 



6EL 



= 2 cr pp c« N (r) 



where x — 



4 E y c 2 
4E 2 + ra 2 c 

7 



qAr,Ejfi) 

"CRpW [bm^c 1 
GeV 

\2 



(17) 



GeV 



(18) 



Owing to Lorentz symmetry, this formula is valid for both lim- 
iting energy regimes, E y » m n o c 2 /2 and E y «: m„o c 2 /2. 
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Because of an incomplete accounting of physical processes at 
the threshold of pion production like the velocity distribution 
of CRp and momentum dependent inelastic CRp-p cross sec- 
tion, eq. d!7t overestimates the number of y-rays for energies 
around E y c 2 /2. 

2.2.3. Dermer's model 

In order to make detailed predictions for the 7r°-decay induced 
y-ray spectrum, more realistic effects near the 7r°-production 
threshold have to be included. This was done by using the 
code COSMOCR originally designed for cosmic ray studies 
by iMiniatil J200ll) . The underlying Av?-isoba ric model wa s 
shown to work well at low proton energies (Stecker 119701) . 
It assumes the CRp-p interaction to be mediated by the exci- 
tation of the A3/2-resonance which subsequently decays into 
two protons and a 7r°-meson. The production spectrum of sec- 
ondary 7r°-mesons is given by a convolution of the normal- 
ized A3/2-isobar mass spectrum represented by a Breit-Wigner 
distribution with the energy dist ribution function. The scal- 
ing model used at high energies ( St ephens & Badhwarl fl98 1 ) 
uses Lorentz invariant cross sections for charged and neutral 
pion production in p-p interactions inferred from accelerator 
data. COSMOCR includes also the contribution of the two 
main kaon decay modes to sec ondary pion spectra (follow- 
ing Moskalenko & Strong 1998) which are K* — > yu* + v^/v^ 
(63.5%) and K* -» n° + (21.2%) where the latter channel 
also contributes to the y-ray source function. 

In order to derive an analytic formula describing the om- 
nidirectional differential y-ray source function over the energy 
range shown in Fig. |2] we keep the behavior of the spectrum 
in the fireball model for E y » rry> c 2 /2 and parameterize the 
detailed physics at the ^-threshold by the shape parameter 5 y 
which smoothly joins the two power laws to the asymptotic ex- 
pansion of the S-function of eq. J17i . yielding 

i «CRd('") 

q y (r, E y ) dE y dV - cr pp c « N (r) t" 1 -p^r < 19 > 

<je v 

3 a y \ GeV / [\ m n o c 1 j \ m n o c z ) 

The scaling behavior in the high-energy limit of Dermer's 
model can be described by a constant pion multiplicity £ = 2 
characterizing the two leading pion jets leaving the interac- 
tion site in direction of the incident protons diametrically and 
carrying the high longitudinal momenta owing to Lorentz con- 
traction of the interacting nuclei in th e center of mass s ystem 
and Heisenberg's uncertainty relation (Nachtmann 1990). This 
assumption of constant pion multip licity of the scaling m odel 
is in contrast to the fireball model ( Mann heim & Schlickeise J 
11994 . which assumes a state of hot quark-gluon plasma in 
thermal equilibrium after the hadronic interaction subsequently 
ablating pions with multiplicities £j0 =s [(E p - £ , t h)/GeV] 1/4 , 
where = 1 .22 GeV denotes the threshold energy for pion 
production. 

The y-ray source function peaks at the energy of nyi c 2 /2 
67.5 MeV. It is well known, that the asymptotic slope of the y- 
ray spectrum, characterized by its spectral index a y , reproduces 
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-0.2 r 



io" 4 io" 2 io° io 2 io 4 

E y [GeV] 

Fig. 2. a) The omnidirectional (i.e. integrated over An solid 
angle) differential y-ray source function q y (E y ) normalized 
by the target number density n^(r) and CRp normalization 
«CRp(r) m order to be independent of the spatial dependence 
of any specific model. The dotted lines show the simulated y- 
ray spectra while the solid curves represent our models given 
by eq. (II 9i with the spectral indices from top to bottom, a v e 
{2.0, 2.4, 2.7, 3.0, 3.5). b) Relative deviation of our analytic 
approach to simulated y-ray spectra. 



the sp ectral index of the population of CRp, a y = a p dPermerl 
1986b). This is again in contrast to the fireball model which 
predicts a steeper asymptotic slope in the y-ray spectrum for 
a p > 2, amounting to a y = 4/3 (a v - 1/2). In the following we 
restrict ourselves to Dermer's model because it is better moti- 
vated by accelerator data. 

By comparing the logarithm of the y-ray source function 
of eq. d!9i to numerically calculated spectra using COSMOCR 
we recognized that the influence of the detailed physics at the 
threshold together with the kaon contribution can be modeled 
in our semi-analytic approach in eq. dl9l by self-consistent 
scaling relations for the shape parameter 6 y and the effective in- 
elastic p-p cross section <x pp including the kaon decay modes. 
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The shape parameter 6 y scales with the spectral index of the 
y-ray spectrum as 



6 y = 0.14a: L6 + 0.44, 



(20) 



which models the functional behavior of the spectrum (com- 
pare Fig.|5Jl. The effective cross section <x pp also depends on a y 
which can be modeled by 



pp 



32 x (0.96 



+ e 



4.4-2.4 



a A mbarn. 



(21) 



On the one hand, the enhanced contribution to the normal- 
ization of the y-ray source function q y {E y ) for fiat spectral 
indices a y is due to the larger contribution of the channel 
p + p — > K* + X relative to p + p — > n ± + X for larger en- 
ergies, approaching asymptotically a value of 27 % at energies 
larger than 1 TeV (Miniati 2001) which we did not account 
for a priori in our simple model. Secondly, this scaling behav- 
ior also includes higher order contributions to the effective pion 
multiplicity for harder spectra characterized by a lower spectral 
index a y . 

The effective description of the spectrum with the smooth 
peak characterized by the shape parameter 5 y starts to fail for 
very steep spectra of a y > 3.5 where relativistic kinematics at 
the threshold plays a crucial role. Then the higher number of 
decaying low energetic 7r°-mesons results in a more concen- 
trated peak on top of the boosted broader distribution of decay- 
ing highly-energetic pions. The lower panel of Fig.[2]shows that 
the relative deviation of the semi-analytic approach of eq. fl!9l > 
to the simulated y-ray spectra amounts below 0.2 for the spec- 
tral range shown in Fig. [2] which is sufficient for the purpose 
of our work. 



2.3. Energy band integrated y-ray luminosity: Analytic 
Ty-Fx scaling relation 

In the following, we derive an analytic T y -Fx scaling relation 
which should serve as an approximate estimate for a given clus- 
ter of galaxies. As a simple scenario we choose the CRp energy 
density to be a constant fraction of the thermal energy density, 
scRp(r) = Xcr p s±(r). However, this is not a fundamental con- 
straint for this scaling relation. Any other spatial dependency 
for the CRp scaling parameter Xcr p may be substituted instead 
of the assumed one. 

The bolometric X-ray emission of the hot thermal intra 
cluster electron s is given by the cooling fu nction for thermal 
bremsstrahlung (Rvbicki & Lightman 1979), 



A x [Bo(r),r e (r)] = A « e (r) 2 jkT e (r) 

In V /2 2 5 7T£> 6 £/, 

, 3 m, 



with Aq = 



3 hm e c 



t3I Z 2 g B (T e ) 



6.62 x 10~ 24 erg s" 1 cm 3 keV 



-1/2 



(22) 
(23) 
(24) 



where n e is the electron number density, T e the temperature, 
oftar is the nucleon density in the ICM relative to the electrons 



for primordial element composition, Z the charge number 1 and 
gs - 1.2 is the frequency and velocity averaged Gaunt factor. 

In order to obtain the integrated y-ray source density A y for 
pion decay induced y-rays the y-ray source function q y {r, E y ) 
in eq. dl9t can be integrated over an energy interval yielding 



A 7 (r,Ei,E 2 ) 

where A y (a p ) 
N y (a p ) 

and Xj 



rE 2 

I dE y q y (r, E y ) 
Je, 

A y (a p ) N y (a p ) X CRp n 2 e (r) kT e (r) , 
cr pp (Qfp) d e 



GeV m n 



[GeWj 



& 



"257' ~2$T j 



(Op - 1) 



1 + 



2E t 



2S„ 



for ie {1,2}. 



Here we introduced the abbreviation 



[/(*)]£ =/(*2)-/(*l). 



(25) 
(26) 
(27) 

(28) 
(29) 

(30) 



Assuming Dermer's model the y-ray spectral index scales as 
a y = a p in contrast to the fireball model where a y = 4/3 (a p - 
1/2). The shape parameter 6 Y is given by the a y -5 y scaling 
relation in eq. J20t which strictly holds for Dermer's model, 
but should also be valid for the extended fireball model. 

Comparing the integrated y-ray source density A Y {E\ , Ei) 
of eq. (I25> to that of thermal bremsstrahlung (eq. (I22» we ob- 
tain an analytic Ty-Fx scaling relation for the ratio of y-ray 
fluxes T y and bolometric X-ray fluxes Fx, 



Ty(E l <E<E 2 ) Ay(a p )N y (a p ) I (kT e ) 



F^ 1 erg- 1 



A keV~ 1/2 erg- 1 



keV 



1/2 



Xcrp, (31) 



where the prefactor is appropriately scaled yielding a dimen- 
sionless number which consists of A y {a p ) (eq. \21Y ). N y (a p ) 
(eq. J28» . and Ao (eq. ( I23t ). The Ty-Fx ratio scales linearly 
with the scaling parameter Xcr p given by eq. @ and is in- 
dependent of the underlying cosmology, however not of red- 
shift due to the K-correction. Inferred values for the expe cted 
y-ray fl ux T y are consistent with those obtained by E nBlin et"all 
( 1997) for the spectral index of our Galaxy a y = 2.7. 

2.4. Stationary spectrum of hadronically originating 
secondary electrons 

This section is bas ed on a formalism developed in 
lDola g& EnBlin (2000). The steady-state CRe spectrum is gov- 
erned by injection of secondaries and cooling processes so that 
it can be described by the continuity equation 



_d_ 

BE, 



(E e (r, E e )f e (r, E e )j = q e (r, E e ) . 



(32) 



1 Setting Z 2 = 1 in eq. I23i is correct for a plasma of primordial 
element composition which consists of hydrogen and helium only, be- 
cause (wnZ 2 ) = «n in this case. This is a reasonable approximation 
owing to the small contamination of heavier elements in the ICM. 
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For E e (r, p) < this equation is solved by 

If 00 , 
\E e (r, E e )\ J Ec 



(33) 



The cooling of the radio emitting CRe is dominated by syn- 
chrotron and inverse Compton losses giving 



- E e (r, E e ) = 



4<T T C 



3 m 



2^4 



8tt 



8tt 



(34) 



where cr x is the Thomson cross section, B(r) is the local mag- 
netic field strength and #c MB /(8 n) is the energy density of the 
cosmic microwave background expressed by an equivalent field 
strength Bcmb = 3.24 (1 + z) 2 pG. The CRe population above a 
GeV is therefore described by a power-law spectrum 



Mr,E e ) 

«CRe(f) 



«CRe(r) / E, 



\ GeV/ 



GeV 

2 7 7rl6 _(Q '^ 1) A mod o-ppm^c 4 «nW "CRpW 



o-jGeV B 2 {r) + B\ 



2 

CMB 



a e = 



a-p + 1 in Dermer's model , 



3(f) 



•(OT.-l) 



in the fireball model , 
in Dermer's model , 
in the fireball model . 



(35) 
(36) 

(37) 

(38) 



For the sake of consistency, we use Dermer's model through- 
out the paper where the effective cross section cr pp is given by 
eq. Mil in contrast to the fireball model where cr pp = 32 mbarn. 
The approach of the scaling relation of eq. J2 1 > is approxi- 
mately valid for CRe although the decay channels of charged 
kaons provide a stronger contribution to the branching ratio 
relative to 7r°-mesons resulting also in slightly higher injection 
rates for electrons and positrons. Differences in normalization 
and radio brightness morphology due to the different models 
governing the CRp-p interaction are small and irrelevant for 
our conclusions. 



2.4.1 . Synchrotron emission of secondary electrons 

The synchrotron emissivity j v at frequency v and per steradian 
of a power law distribution of CRe (eq. (I35» in an isotropic 
distribution of magneti c fields and electrons w ithin the halo 
volume (eq. (6.36) in lRvbicki & Lightman ll979l) . is obtained 
after averaging over an isotropic distribution of electron pitch 
angles yielding 



(39) 



Mr) = c 2 (a e )hc Re (r)B(rr^\- 
with d = 3 e GeV 2 /(2 nml c 5 ) , 



y/3n 



a e + 



r(¥)r(¥)r(^) 



32 n m,c 2 a e + 1 



r(fl 



(40) 



where r(a) denotes the T-function dAbramowitz & St egun 
1 19651) and a v = (a e - Fj/2 = a p /2 in Dermer's model. In our 
models the magnetic field B(r) was assumed to be spherically 



symmetric on cluster core scales and a nd to follow the electron 
density n e (r) fctolag et alll999ll200ll) : 



B(r) = B 



n e (r) 



«e(0) 



(41) 



where Z?o an d a B are f ree parameters in our model. Assuming 
the radio emissivity j v (r) in eq. fl39l > to be only a function of ra- 
dius, then the line of sight integration yields the surface bright- 
ness of the radio halo 



Sv(rx) = 2 



r R Mr)r 



(42) 



2.4.2. Inverse Compton emission of secondary 
electrons 

The source function gic owing to IC scattering of cosmic mi- 
crowave background (CMB) photons off an isotropic power 
law distribution of hadronically originating CRe (eq. \35Y ) is 
(derived from eq. (7.31) in[Rvbicki & Li ghtmaij l 19791 in the 
case of Thomson scattering), 



qic(r, E y ) 
/ic(«e) 

and q{r) 



HK >HCK ejy^J [kTcMBj 

2 a " +3 (a 2 e + 4 a e + 11) 
(a e +3) 2 (a e + 5)(a e + 1) 
„ / a e + 5 \ la e +5 

xr (^)n^ 

8 n 2 r\ h C R e (r) (JcTcmb) 2 
h 3 c 2 



-(<*.+« 



(43) 
(44) 
(45) 
(46) 



where a v = (a e - l)/2 denotes the spectral index, r e = 
e 2 /(m a c 2 ) the classical electron radius, ((a) the Riemann 
function dAbramowitz & Stegunlll965l) . and ncReW is given 
by eq. d36i . After integrating over the IC emitting volume in 
the cluster we obtain the particle flux T{E y ) (see eqs. l|4} and 
d5}). The same CRe population seen in the radio band via syn- 
chrotron emission scatter CMB photons into the hard X-ray 
regime. In the y-ray spectrum, there is a point of equal con- 
tribution of the IC spectrum of the CRe showing a decreasing 
slope of -a v — 1 = — Qfp/2 - 1 (assuming Dermer's model) and 
the pion decay induced y-ray spectrum being characterized by 
the rising slope a y - a p (see eq. dl9». In the high energy limit 
(E y » m^o c 2 /2), the pion decay induced y-ray spectrum de- 
clines with a slope of -a y = —a p which is the same as the IC 
emission for a p — 2 and slightly steeper for larger values of a p 
(for illustration, see Fig.|3J. 

2.5. Summary and outline 

In the previous sections we developed an analytic description to 
compute the neutral pion decay induced y-ray spectrum from a 
CRp population over a broad range of y-ray energies extending 
from below MeV up to TeV. Moreover, we presented a for- 
malism of calculating the synchrotron and IC emission from 
a stationary population of CRe resulting from hadronic CRp 
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Fig. 3. The simulated differential flux of y-rays from Perseus 
reaching the Earth. Shown are upper limits of the IC emission 
of secondary CRe (power-laws, assuming zero magnetic field) 
as well as pion decay induced y-ray emission (represented by 
broad distribution centered on £p ea k — 67.5 MeV). The nor- 
malization of the spectra differing in their values of the CRp 
spectral index a y - a p (Dermer's model) depends on the as- 
sumed scaling between CRp and thermal energy density. We 
fix this scaling parameter Xqr p assuming the isobaric model 
(Sect. lXXTl by comparing the inte grated flux above 100 MeV 
to EGRET upper limits (see Reim eret all2003l) . 

interactions of the thermal gas of the ICM. Assuming a con- 
stant scaling between kinetic CRp energy density and thermal 
energy density of the ICM we derived an analytic Ty—Fx scal- 
ing relation for the ratio of y-ray flux to bolometric X-ray flux 
to obtain observationally promising cluster candidates for con- 
straining the CRp population. In order to obtain reliable flux 
estimates we are going to introduce in the following three pos- 
sible spatial distributions of the CRp, whose population is ei- 
ther in fractional pressure equilibrium with the thermal parti- 
cle population (as assumed for selecting the clusters), experi- 
enced adiabatic compression during the formation of the cool- 
ing flow cluster or is shaped by diffusion away from a central 
source of CRp. By modeling the y-ray emission of these par- 
ticular clusters and comparing to EGRET upper limits we are 
going to present bounds on the CRp population. Furthermore, 
we will derive upper bounds on the CRp population by radio 
synchrotron emission of hadronically originating CRe and will 
compare azimuthally averaged radio brightness profiles of the 
the Perseus radio-mini halo and the radio halo of Coma. 

3. Astrophysical application to nearby clusters of 
galaxies using multi-frequency observations 

3.1. The expected spectral index a p 

The spectral index of the ICM CRp population a p is not well 
constrained by observations. However, becau se galaxy clusters 
are able to store CRp for cosmological times JVolket all 1993 
Enfi lin et aDll997t iBerezinskv et aljfl997l) the spectral index 



of the global CRp population (allowing for spatial differentia- 
tion) is expected to be that of the injection process, if no re- 
acceleration processes modified the spectrum after injection. 
We discuss briefly different possible CRp sources in galaxy 
clusters: 

Structure formation shock waves have generated most of the 
thermal energy content of galaxy clusters. Therefore, it is plau- 
sible to assume that they also produced most of the CR en- 
ergy of clusters. Shock acceleration is able to produce momen- 
tum power-law particle distributions characterized by a spectral 
index, which is in the test-particle picture of non-relativistic 
shock acceleration 



_ R + 2 
a ' ms ~ ~R~T ' 



where R < 4 is the shock compression factor. The lowest spec- 
tral indices are therefore generated by the strongest shocks, 
whic h are preferentially found in peripheral regions of the clus- 
ters dOuilis et alJll998tlMiniati et aljboOOl) . Thus, harder CRp 
populations (a^ = 2.0. ..2.5) are mostly injected into the out- 
skirts of clusters. However, motion of the ICM gas transports 
them efficiently into the cluster centers jMiniati et all2 001b). 
Injection by radio galaxies: Active galactic nuclei (AGN) 
are able to produce large amounts of relativistic plasma. The 
composition of this plasma is not known, however, the pres- 
ence of CRp is assumed in and supported by many papers. 
The energetics of AGN is sufficien t to inject a signific ant CRp 
population into the thermal ICM jEnBlin et all ll 9971 [T9 98b; 
IColafrancesco & Blasilfl998l lBlasilfl999HWu et alJl2000l) pro- 
vided CRp are present in the radio plasma and are able to leave 
it on cosmological short timescales. If the radio plasma releases 
all its CRp, a moderately flat injection spectrum can be ex- 
pected (say ffj n j w 2.5) since radio emission from radio galaxies 
indicates flat CRe spectra. If, however, only a small fraction of 
the CRp is able to leave the radio plasma diffusively, an even 
flatter spectrum (say a u ^ « 2.2) can be e xpected due to increas- 
ing escape probability with momentum (EnBlin 2003). 
Supernova Remnants (SNR) are known to be able to produce 
flat (ai n j » 2.4) CR populations and they are believe d to be the 
main CRp source of our galaxy JSchlickeiserl2002l and refer- 
ences therein). The reason for the steeper (a p m 2.7) galactic 
CRp spectrum is thought to be the momentum dependent es- 
cape probability from our Galaxy. Thus, the spectrum of CRp 
escaping from galaxies should be flat again (a^ * 2.4). The 
spectrum injected into galaxy clusters could be even flatter, if 
termination shock waves of the galac tic w i nds ar e able to re- 
accelerate them, as proposed by Volk et al. ( 1996). 

The CRp population in galaxy clusters which is able to 
interact with the thermal gas and thus to produce observable 
signatures will be a mixture of contributions of the different 
sources, modified by acceleration and energy loss processes. 
In order not to rely too much on a specific physical picture, 
we discuss simplified models, which should be able to capture 
many typical situations. 
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3.2. Spatial distribution of cosmic ray protons in 
cooling flow clusters 

In the following, we introduce three models for the spatial dis- 
tribution of CRp within clusters of galaxies. The origin of the 
CRp population is not specified in the first two models, but 
the CRp may be accelerated by shock waves of cluster merg- 
ers, accretion shocks, or result from supernova driven galactic 
winds. In contrast to that we explore in the third model the 
diffusion process of CRp away from a central AGN. Since it 
is unclear how CRp are distributed spatially in detail, we in- 
vestigate here three different scenarios which should serve as 
toy models. We pursue the philosophy of estimating physical 
parameters from observationally obtained electron density and 
temperature profiles by using simplified model assumptions for 
the CRp population. In doing so we try to minimize the di- 
mensionality of parameter space as much as possible in order 
to track the main physical processes by means of analytically 
feasible methods and not to rely upon too many assumptions. 
Therefore the presented CRp profiles which are based on the 
assumption of spherical symmetry should not be interpreted as 
a precise estimate of the CRp population but rather as a plausi- 
ble spherically averaged scenario. 



3.2.1 . Isobaric model of CRp 

In this model we assume that the average kinetic CRp energy 
density £cr p (/) is a constant fraction of the thermal energy den- 
sity £ t h0) of the ICM 



ScRpW = XcRp e th (r) . 



(48) 



This distribution might be maintained even in the case of a 
cooling flow cluster by mixing and ongoing turbulent CRp dif- 
fusion processes exerted by relativistic plasma bubbles rising in 
the gravitational poten tial of the cluster due to buoyant forces 
JChurazov et al.l200l[ and references therein) which possibly 
leads to fractional pressure equilibrium with the thermal parti- 
cle population. 

3.2.2. Adiabatic compression of CRp 

Here we assume the CRp population to be originally isobaric to 
the thermal population but to become adiabatically compressed 
during the formation of the cooling flow while it did not relax 
afterwards. The phase space volume stays constant during this 
transition and the momenta and volumes scale according to 



Pp 



1/3 



Pp 



c l/3 



Pv. 



(49) 
(50) 



Here the compression factor C = C(r) = (n' e /n e )(r) has been 
introduced, which is larger than unity within cooling flows. 
Provided that the electrons have been in hydrostatic equilib- 
rium during this transition, this implies C{r) = T c \ ustel /T^,{r), 
where r c i uster denotes the electron temperature in the outer core 
region. This transformation implicitly assumes that the ratio of 



the CRp number densities before and after the adiabatic com- 
pression equals that of the electron population. If the differen- 
tial number density distribution of the CRp population may be 
described by a power-law in momentum p p , then after adiabatic 
compression of CRp the functional shape of their distribution 
remains unchanged, however shifted according to 



/»') 



l CRp 



(r')c 



(51) 
(52) 



GeV \GeVJ 

"CRp^O = «CR P [r'(r)] C(r') ( ^ +2)/3 

The normalization «cRp(f) is chosen in such a way that the 
kinetic CRp energy density makes up a constant fraction of the 
thermal energy density prior to cooling flow formation and is 
described by a scaling parameter Xcr p 

fifcupC) = X CRp £thW -» ecRpC*-') = x CR P ( r ') fi th(r') . (53) 

After adiabatic compression of CRp this scaling parameter has 
thus changed to 

X' CRp (r) = C (a > +2)/ \r)X CRp . (54) 

Since any hadronically induced emissivity scales with X' CRp we 
obtain the following relation, 



c( r ) _ c (a p +2)/3^ -isobaric^ 



3.2.3. Diffusion of CRp away from a central AGN 



(55) 



Many galaxy clusters - especially those with a cooling flow 
- harbor a central galaxy, which often exhibits nuclear activity. 
The relativistic plasma bubbles produced by the AGN may con- 
tain relativistic protons, which can partly escape into the ther- 
mal ICM. Most of the CRp that have been injected into the clus- 
ter center are either diffusively transporte d into the surround- 
ing IC M (as assumed by IColafrancesco & Blasil |l9 98 ; Blasi 
1 999) or form relativistic bubbles which rise in the gravitational 
potential of the cluster due to buoyant forces (Churazov et al. 
2001, and references therein). An argument in favor of a sig- 
nificant central CRp injection into the ICM is the much more 
efficient escape of CRp from the magnetic confinement of the 
radio plasma bubble during the very early stages due to the bub- 
bles higher geometrical compactness and and expected stronger 
turbulence level (EnBlin 2003). In addition to this, any galactic 
wind from a central galaxy will also inject CRp into the cluster 
center. In order to treat these diffusion processes analytically 
one has to distinguish between clusters containing a cooling 
flow region or not. In the first case CRp diffusion will shape 
their emission profiles owing to the peaked cooling flow pro- 
files while the emission strength in non-cooling flow clusters is 
mainly governed by the effective injection timescale. 

Cooling flow clusters: The transport of CRp through the 
ICM is diffusive, with a diffusion coefficient /c(r, p) which in 
general may depend on momentum and position. For illustra- 
tion we use 



K(r,p) = Ko(r 



I pc p 
; \GeV/ ' 



(56) 



with Kq ~ 10 cm s being plausible values. By using 
this ansatz, we ignore likely deviations of the diffusion coef- 
ficient from eq. ( I56> in the mildly relativistic regime because 
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these CRp are also not constrained by observations of their 
hadronic interactions. The coefficient a^s describes the mo- 
mentum dependence of the diffusion and is expected to be 
<*diff ~ 5 for active CRp diffusion in a Kolmogorov-like small- 
scale magnetic turbulence spectrum and a&g ~ for pas- 
sive advective transport in a turbulent flow. In the latter case 
«(p) = Vturb^turb/3 ~ 10 29 cm 2 s _1 , where Vtmb ~ lOOkm/s and 
^turb ~ 10 kpc are the turbulent velocity and coherence length, 
respectively. 

In a stationary situation, which is a valid approximation for 
timescales longer than the typical CRp diffusion timescale in 
the case of a stationary or short-term intermittent CRp source, 
the CRp distribution functions is given by 



fp(r,P P ) 
where Q v (p) 



dr' 



Qp(p) 



Q(P) = 

4?r Jo K(r',p)r' z 4nK(p)r 
Qp.qc ( PpC 
GeV 



I GeV/ 



(57) 
(58) 



is the averaged CRp injection rate of the central source. We 
assume it to be a power-law in momentum with spectral in- 
dex ainj , which in general is not identical to the spectral index 
of the CRp population within radio plasma si nce the escap e 
fraction is expected to depend on momentum (EnBlin 2003). 
For the last step in eq. (I57t we assume for simplicity the dif- 
fusion coefficient to be independent of position. Possible mod- 
els of spatial distributions for the diffusion coefficient depend 
strongly on many unknown quantities such as the dominant 
diffusion mechanism (active diffusion versus passive advective 
transport), the velocity field, the turbulence scale, and the topol- 
ogy of the magnetic field, only to mention a few. Therefore 
we are unable to guess a realistic profile for Ko(r) without en- 
larging the accessible parameter space tremendously. However 
since we expect the diffusion coefficient not to change dramati- 
cally over the cooling flow scale and since the distribution func- 
tion / p (r, p p ) is sufficiently steep in radius (eq. i57\ ) our results 
should be approximately correct. The total CRp luminosity of 
the source can be estimated from eq. i57\ to be 



In order to obtain a realistic estimate for the diffusion vol- 
ume to be considered, the relevant length scale needs to be 
taken into account. We define the characteristic scale R^g by 
calculating the second moment of the time-dependent distribu- 
tion function of the first particles released by the source, yield- 
ing 



/ -1/2 / P c V 

fldiff = y2n dim f inj K(p) a 80/z 70 ' kpc [^-^ J 



ffdilT/2 



(61) 



where n& m = 3 denotes the number of spatial dimensions. 
Here we assume a typical lifetime of fj n j = 3 lui Gyr and 
a:o — 10 29 cm 2 s _1 . Beyond this scale R^itf there can be a CRp 
population resulting from diffusion away from the central AGN 
which is however exponentially suppressed. Because the y-ray 
luminosity resulting from hadronic CRp interactions scales as 



i; r oc47r J drr 2 n C R V (r)n e (r) cc j dr i 



(62) 



we always obtain centrally peaked y-ray profiles, since the 
cooling radius is smaller than the diffusion scale, r C] < R^is, 
and p > 1/3 within cooling flow regions (compare Table 
Thus, the y-ray luminosity is only weakly dependent on R^ as 
long as it reflects the correct order of magnitude. 

In this work we constrain «cRp,o with the aid of y-ray ob- 
servations of galaxy clusters. From these constraints limits on 
the averaged CRp luminosity escaping from the radio plasma 
of the central galaxy can be derived using 



L CRp 4 n ot p c 2 «cr p ,o fc 7 o kpc / m p c 



70 

2(ainj - 1) 
ai„j - 2 3 - a. 



GeV 



inj 



(63) 



where we again ignored any possible low-energy spectral cut- 
off, since it can be included a posteriori with the help of Fig.^ 
As a rough estimate we find numerically 



^CRp = 



2 r\ 1 2 \ l _Q, in 



2(otoj-l) \GeV 



^inj 2 3 Q'inj 



(59) 



Within our model, the CRp distribution function within the 
thermal ICM can be written as 



fp(r,p v ) = 



«CRp,0 c 

GeV 



'ho 



kpc J 



(Pp c Y 
\ GeV/ 



(60) 



where a p = orj n j + a&ff ■ 



2 This seems to be in contradiction to the identity of injection and 
equilibrium spectral index for a system without escape claimed in 
Sect. 13. II Formally, we had to include particle escape from the galaxy 
cluster in order to be able to have a finite steady state solution of the 
diffusion problem, as given by eq. <57> . In the realistic case of a finite 
age of the system the stationary solution is only approximately valid 
in the center of the galaxy cluster. However, only there exists a suf- 
ficiently high target density to detect the CRp population. Therefore, 
although we use a poor description of the CRp profile on large-scales, 
the estimated y-ray fluxes should be sufficiently accurate. 



^CRp 



L(a inj ) 10 43 ^ 1/2 ergs- 1 



( «! ) 

U0 29 cm 2 s-V 



«CRp,0 



10- 6 h^cm- 



(64) 



with L(a; n j) = 6.1, 2.2, 1.6, and 1.7 for a m j =2.1, 2.3, 2.5, and 
2.7, respectively. In Sect l3.3.4l we analyze these constraints for 
our cluster sample in more detail. 

Non-cooling flow clusters: In transforming the above con- 
siderations on diffusion length scales to the case of non-cooling 
flow clusters we point out the following differences: In non- 
cooling flow clusters the core radius is normally larger than 
the diffusion scale, r c > Rdm, over which the electron density 
varies only slightly. Thus, a stationary solution to the diffusion 
equation is not applicable in the case of a flat target profile. It 
follows that the volume integrated y-ray spectrum does not de- 
pend on the diffusion coefficient but only on the injection time 
fjnj of CRp into the ICM of the cluster core. We therefore adopt 
a modification to the diffusion model for non-cooling flow clus- 
ters. The averaged CRp luminosity of the central galaxy reads 
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in this context 

N C Rp 



^CRp - 



2 (a, 



inj 

2 3 



1) \ GeV 



(65) 



Here Ncr p denotes the integrated number of CRp being in- 
jected into the 1CM of the cluster and a m j = a p , because there is 
no diffusion induced spectral steepening simply due to the fact 
that the even more energetic CRp which are still significantly 
contributing to the y-ray flux in the EGRET energy band are 
not able to leave the central core region within a reasonable 
timescale. 

3.3. Constraining the population of CRp by the 
integrated flux ofy-rays in different clusters 

3.3.1 . Cluster sample 

Applying the T y -Fx scaling relation (eq. (13 \\ ) and taking 
bolometric X-ray fluxes from bavid etalJ (ll993) while fixing 
^CRp = 0.01 Xo.oi and a p = 2.3 we estimated y-ray fluxes 
7>,est(> 100 MeV) for the spectral sensitivity of EGRET in or- 
der to choose our cluster sample (see Table [2}- Inferred val- 
ues for the estimated y-ray flux Ty^st by means of the T y - 
Fx scaling relation sensitively depend on the bolometric X- 
ray luminosity of the particular cluster such that values for 
^y.est in Table |2 represent a rough estimate. A detailed mod- 
eling using density and temperature profiles will be described 
later on in Sect. 13.3.21 in order to obtain upper limits on the 
CRp population. By comparing y-ray fluxes T y obtained from 
these two different methods we recognized an inconsistency 
for the Virgo and Centaurus cluster: This discrepancy is ex- 
plained b y a too sma l l apert ure of the X-ray experiments ana- 
lyzed bv lDavid etaD i 19931) giving rise to an underestimation 
of the X-ray flux of these two nearest clusters in our sample 
(zvirgo = 0.0036 and zcentaurus = 0.0114) and therefore an un- 
derestimate of :F y ,est for these two clusters. Moreover, we no- 
ticed a systematic discrepancy of the order of 50% between the 
different methods in cooling flow clusters which is due to an 
insufficient accounting for the radial temperature variation in 
eq. (£T}. 

Parameters of electron density profiles n e {r) of our cluster 
sample are given in Tablen wnere the clusters are ordered ac- 
cording to their property of containing a cooling flow (upper 
part) or not (lower part). Note that the parameters are subject 
to different formulae d66i and {671 . 



§44) 



n e (r) = 



1=1 

A[Te(0)] 



1/2 



(66) 



(67) 



The last equation i67\ follows from deprojection of X-ray 
surface brightness profiles which are represented by double 
P models. The derivation of this deprojection is given in 



AppendixlAl For simplicity and consis tency with th e X-ray sur- 
face brightness profiles given in lMohr et alJ i 19991) we ignored 
the weak dependency on T e (r) in eq. i67i . 

In order to model the temperature profiles T e {r) for our 
cooling flow cluster sample we applied the universal tempera- 
ture profile for relaxed clusters proposed bv lAllen et alJ d200lh 
to data taken from the literature, 



T e {r) = T Q + (T\ - T ) 



1 + 



' temp 



(68) 



This equation matches the temperature profile well up to radii 
of ~ 0.3 r v j r , which is sufficient for our purposes since we are 
especially interested in the core region of clusters. The param- 
eters of the temperature profile for particular cluster are given 
in Table 

3.3.2. Simulated y-ray flux normalized by EGRET 
limits: The case of Perseus cluster 

The volume integrated omnidirectional differential y-ray 
source function Q Y (E y ) can be obtained by integrating eq. dl9l . 
We integrated the volume out to a radius of 3 hzl Mpc which 
corresponds to the characteristic distance where the simple 
j6-model of electron densities breaks down due to accretion 
shocks in clusters. The integration kernel q 7 (r, E y ) scales lin- 
early with «cRp(f) (as shown in eq. j!9H which is obtained by 
solving eqs. @ and dlOi . By comparing the integrated y-ray 
flux above 100 MeV, TJ > 100 MeV), to EGRET upper limits 
("see lReimer et all2 003). we constrain the CRp scaling param- 
eter XcRp- The inferred value for Xcr p in the Perseus cluster 
normalizes the differential y-ray flux 



dT Y = Q y (E y ) 
dE v ~ AttD 2 



(69) 



in Fig. The 7r°-meson decay induced distinct spectral sig- 
nature resulting in the peak at a y-ray energy of m n o c 2 /2 ^ 
67.5 MeV can be clearly seen. 

Figure [3] shows also upper limits on the differential y-ray 
flux owing to IC emission of hadronically originating CRe rep- 
resented by power-laws. The IC spectra are computed by means 
of eq. J43l > for different spectral indices a p and zero magnetic 
field. Non-zero magnetic fields can be included since the IC 
spectra scale according to B^, MB /(B(r) 2 + #cmb) ( see e 1- OH 1 ) 
which results in a lower normalization. 



3.3.3. Results on the scaling parameter X C r p using 
y-ray observations in different clusters 



By employing the technique described in Sect. 13.3.21 we con- 
strained the CRp scaling para meter Xcxr, using E GRET up- 
per limits of the y-ray flux by Rei mer et alJ ( 120031) . As de- 
scribed in that section, we infer the y-ray flux of this clus- 
ters originating from within a sphere of radius 3 hzl Mpc. 
Owing to the vicinity of the Virgo cluster this maximum ra- 
dius subtends an angle on the sky which is larger than the 
width of the point spread functio n of the EGRET in strument 
(6/ max = 5.8° [£,7(100 MeVir 534 . lReimer et alJll2003l) l. Thus 
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Table 1. Parameters of electron density profiles n e {r) of our cluster sample (central densities «, are subject to different formulae 
d66t and (I67» . The cluster are ordered according to their property of being a cooling flow cluster (upper part) or a non-cooling 
flow cluster (lower part). 



Cluster 


z 


»i 

[hjQ cm- 3 ] 


[Kjl kpc] 




112 

[h% cm- 3 ] 


Uho kpc] 




Equation 


References 


A85 


0.0551 


3.08 x 10 2 


45 


0.662 


3.87 x 10~ 3 


226 


0.662 




(a), (b) 


A426 ("Perseus') 


0.0179 


4.6 x 10~ 2 


57 


1.2 


4.79 x 10~ 3 


200 


0.58 


l66l 


<c\ (eft 


A2199 


0.0302 


3.37 x 10~ 2 


29 


0.663 


7.17 x 10" 3 


116 


0.663 


|67} 


(a), (b) 


A3526 (Centaurus) 


0.0114 


8.05 x IQ 2 


8.6 


0.569 


3.65 x 10- 3 


99 


0.569 


03 


(a), (d) 


Ophiuchus 


0.0280 


1.71 x 10~ 2 


56 


0.705 


7.47 x 10~ 3 


190 


0.705 


|67} 


(a), (e) 


Triangulum Australis 


0.0510 


7.31 x 10~ 3 


151 


0.816 


2.63 x 10~ 3 


444 


0.816 


|67} 


(a), (f) 


Virgo 


0.0036 


1.5 x 10-' 


1.6 


0.42 


1.3 x 10~ 2 


20 


0.47 


167} 


(g), (h) 


A 1656 (Coma) 

A2256 

A2319 

A3571 


0.0231 
0.0581 
0.0557 
0.0391 


3.4 x 10~ 3 
3.57 x 10" 3 
7.35 x 10~ 3 
9.37 x 10~ 3 


294 
347 
152 
124 


0.75 
0.828 
0.536 
0.61 








|66} 
166} 
166} 
166} 


(i), (d) 
(a), (d) 
(a), (d) 
(a), (d) 



(a) iMohr et al 



McHa rdv et al. 



Il999h . ( b) lOeeerle & Hill l200ll) . (c) Ichurazov et alJ l2003h. (d) Istruble & Rood fl99' 
<198ll) . (g) lMatsushita et alJl2002l) . (h) lEbeling et alJ<1998h . (i) lBriel et alHl992l) 



(e) lLahav et alJ ll989h . (f) 



Table 2. Parameters of temperature profiles T e (r) of our cluster sample. The estimated y-ray flux Ty.^v (> 100 MeV) was calcu- 
lated using the T y -Fx scaling relation (eq. (13 IV ) with a p = 2.3 and bolometric X-ray fluxes from lDavid et al.l jl993h . Note that 
T y scales linearly with Xcr p which was set to Xcr p = 0.01 Xo.oi in this table. The range for T y ^\ reflects the temperature spread 
in cooling flow clusters between the central temperature To and the peripheral temperature T\ . 



kT m r tcmp 7>,est (> 100 MeV) 

Cluster Experiment [keV] [keV] [/Ci kpc] rj [X 00i 10 -10 cm -2 s _1 ] References 

A85 BeppoSAXMECS 5.5 9.0 312 2 1.8 ... 2.3 (a), (b) 

A426 (Perseus) .... XMM-Newton MOS 3.0 7.0 94 3 19.1 ...29.2 (c) 

A2199 Chandra ACIS 1.6 4.3 21.5 1.8 1.3 ... 2.2 (d), (e) 

A3526 (Centaurus) ASCA GIS 2.2 4.0 22 3 2.2 .. . 2.9 (f) 

Ophiuchus ASCA GIS 12.8 22.0 (f) 

Triangulum Australis ASCA GIS 10.3 4.8 (f) 

Virgo XMM-Newton PN/MOS 1.0 3.0 13.5 1 3.2 .. . 5.6 (g) 

A1656 (Coma) XMM-Newton MOS 8.3 13.1 (h) 

A2256 Chandra ACIS 6.7 2.3 (i) 

A2319 ASCA GIS 9.7 5.4 (f) 

A3571 ASCA GIS 7.2 4.5 (f) 



(a) llrwin & Bregmanl l200d). (b) Lima Neto et alJ 1200 ll). (c) Ichurazov et alJ 120031) (d) IVoiet et alJ 120021) . (e) Ijohnstone et alJ 120021) . (f) 
IWhitdfcOOrj ). (g) lMatsushita et alJ <2002\ (h) Arnau d et alj<200ll) . (i) ISun et aljj2002^ 



in the case of Virgo we use this smaller integration volume. 
Table [3] shows constraints for Xqr p using the isobaric and the 
adiabatic model of CRp described in Sect. 13.21 Because in the 
adiabatic model the CRp scaling parameter Xcr p is a function 
of radius, the value X^ bMic given in Table prefers to the un- 
primed quantity in eq. 1541 which reflects the outer core re- 
gion of the cluster. For clusters like Perseus (A 426), Virgo, 
Ophiuchus, and Coma (A 1656) we can obtain quite tight con- 
straints on the population of CRp. 

3.3.4. Results on L C r p in the AGN-diffusion model 

The procedure of inferring constraints on CRp diffusing away 
from a central source is mostly sensitive to the CRp population 



of the central cooling flow region rather than the shock region 
in the outer parts of the cluster. In order to constrain the CRp 
density parameter ncR P; o and averaged CRp luminosity Lcr p of 
the central active galaxy in our AGN-diffusion model of cool- 
ing flow clusters we have to calculate the volume integrated 
omnidirectional differential y-ray source function Q y (E y ) (see 
eq. 13). The integration kernel q y {E y ) is proportional to «cRp(f) 
(eq. M9V ) which is obtained by solving eqs. (|6} and d60> . By 
comparing the integrate d y-ray flux ab ove 100 MeV to EGRET 
upper limits (see IReimer et alJ l2003). we constrain the CRp 
density parameter «cr p ,o- In the case of non-cooling flow clus- 
ters we constrain the averaged CRp luminosity Lcr p with the 
aid of the integrated CRp number parameter Vcr p , yielding an 
indirect measure of a combination of the CRp escape fraction 
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Table 3. Upper limits on the CRp scaling parameter Xqr p by comparing the integrated flux above 100 MeV to EGRET upper 
limits assuming a y-ray spectral index in Dermer's model a y = a p . The spatial distribution of CRp is given by the isobaric and 
the adiabatic model of CRp, respectively (see Sects. l3~2~Tl and l3~2. 21 



T y (> 100 MeV) X\*-£- • X[ 



Cluster [10 8 cm 2 s '] a p = 2.1 a p = 2.3 a p = 2.5 a p = 2.7 a p = 2.1 a p = 2.3 a p = 2.5 a p = 2.7 

A85 < 6.32 3.53 1.97 2.09 3.11 2.58 1.41 1.48 2.16 

A426 (Perseus) .... < 3.72 0.14 0.08 0.08 0.13 0.12 0.06 0.07 0.10 

A2199 < 9.27 6.14 3.42 3.64 5.42 5.74 3.18 3.38 5.00 

A3526 (Centaurus) < 5.31 1.54 0.86 0.91 1.36 1.45 0.80 0.85 1.26 

Ophiuchus < 5.00 0.30 0.17 0.18 0.26 

Triangulum Australis < 8.13 1.93 1.07 1.14 1.70 

Virgo < 2.18 0.18 0.10 0.11 0.16 0.16 0.09 0.09 0.14 

A1656 (Coma) < 3.81 0.45 0.25 0.27 0.40 

A2256 < 4.28 3.15 1.75 1.87 2.78 

A2319 < 3.79 0.86 0.48 0.51 0.76 

A3571 < 6.34 1.85 1.03 1.09 1.63 



Table 4. Cooling flow clusters: Upper limits on the CRp density parameter «cRp,o and average CRp luminosity LcRp of the 
central active galaxy by comparing the integrated flux above 100 MeV to EGRET upper limits assuming a y-ray spectral index 
in Dermer's model a y — a p . The spatial distribution of CRp is is calculated according to the diffusion model of CRp away from 
a central AGN assuming a p — a; n j + a^, where aaiff = 1/3- Non-cooling flow clusters: Upper limits on the CRp number 
parameter Ncr p and average CRp luminosity LcRp without any diffusion induced spectral steepening, i.e. a p = armj. Note that 
LcRp scales in the case of cooling flow clusters with the diffusion coefficient a- while it only depends on the CRp injection time 
finj for non-cooling flow clusters (see Sect. 13.2.31 . 



"cRp.o [h l 7 L 2 cm 3 ] 



L CRp 



\h-l 12 erg s- 1 



CF Cluster 


a p 




2.4 


a P 




2.5 


«p 




2.7 


ffp 




2.9 


a v 




2.4 


a p 






2.5 


ffp = 


2.7 


a p 




2.9 


A85 


6.0 


X 


io- 5 


7.2 


X 


IO" 5 


9.9 


X 


io- 5 


1.3 


X 


10" 4 


5.5 


X 


10 45 


2.7 


X 


10 45 


1.9 x 


10 45 


2.0 


X 


10 45 


A426 (Perseus) .... 


2.4 


X 


IO" 6 


2.9 


X 


io- 6 


3.9 


X 


IO" 6 


5.2 


X 


io- 6 


2.2 


X 


10 44 


1.1 


X 


10 44 


7.4 X 


10 43 


8.1 


X 


10 43 


A2199 


3.2 


X 


IO" 5 


3.8 


X 


io- 5 


5.3 


X 


10 5 


7.0 


X 


io- 5 


3.0 


X 


10 45 


1.4 


X 


10 45 


9.9 x 


10 44 


1.1 


X 


10 45 


A3526 (Centaurus) 


4.3 


X 


10~ 6 


5.1 


X 


10~ 6 


7.1 


X 


10~ 6 


9.4 


X 


10~ 6 


3.9 


X 


10 44 


1.9 


X 


10 44 


1.3 x 


10 44 


1.5 


X 


10 44 


Ophiuchus 


1.5 


X 


io- 5 


1.8 


X 


io- 5 


2.5 


X 


10~ 5 


3.3 


X 


io- 5 


1.4 


X 


10 45 


6.6 


X 


10 44 


4.6 x 


10 44 


5.1 


X 


10 44 


Triangulum Australis 


1.4 


X 


10~ 4 


1.7 


X 


10~ 4 


2.3 


X 


io~ 4 


3.1 


X 


io~ 4 


1.3 


X 


10 46 


6.2 


X 


10 45 


4.3 x 


10 45 


4.8 


X 


10 45 


Virgo 


2.5 


X 


IO" 7 


3.1 


X 


io- 7 


4.2 


X 


IO" 7 


5.6 


X 


io- 7 


2.4 


X 


10 43 


1.1 


X 


10 43 


7.9 x 


10 42 


8.7 


X 


10 42 














NcRp 


















^CRp 




/ 


3/2 -1 / 
70 er 8 S \3 


3^) 








NCF Cluster 


a P 




2.1 


a p 




2.3 






2.5 


ffp 




2.7 


a p 




2.1 


ffp 






2.3 


a„ = 


2.5 


a p 




2.7 


A 1656 (Coma) 


1.6 


X 


10 64 


2.5 


X 


10 64 


3.7 


X 


10 64 


5.1 


X 


1Q 64 


2.7 


X 


10 45 


1.5 


X 


10 45 


1.6 x 


10 45 


2.4 


X 


10 45 


A2256 


1.1 


X 


10 65 


1.8 


X 


10 65 


2.6 


X 


10 65 


3.6 


X 


10 65 


1.9 


X 


10 46 


1.1 


X 


1Q 46 


1.1 X 


10 46 


1.7 


X 


10 46 


A2319 


4.5 


X 


10 64 


7.0 


X 


10 64 


1.0 


X 


10 65 


1.4 


X 


10 s5 


7.5 


X 


10 45 


4.2 


X 


10 45 


4.4 x 


10 45 


6.6 


X 


10 45 


A3571 


2.8 


X 


10 64 


4.4 


X 


10 64 


6.5 


X 


10 64 


9.0 


X 


10 64 


4.7 


X 


10 45 


2.6 


X 


10 45 


2.8 x 


10 45 


4.2 


X 


10 45 



from the radio plasma of the central galaxy and the averaged 
CRp luminosity of this source. 

Upper limits on the CRp density parameter «cr p ,o, number 
parameter of CRp Ncr_ p , and averaged CRp luminosity LcRp of 
the central active galaxy (by means of eq. ( I63» are presented in 
Tabled This shows that within this conceptually simple model 
we are able to put constraints on the averaged CRp luminosity 
LcRp- The limits which are strongest in the case of M87 in the 
Virgo cluster represent conservative bounds since we choose 
the active CRp diffusion scenario resulting in spectral steepen- 
ing of the CRp population. We obtain even tighter limits when 



assuming a passive advective transport of the CRp in a turbu- 
lent flow in which case we infer 

L CRp = L(a p ) 10 42 h-f erg s" 1 ( ^^-t ) (™) 

with L(a p ) = 4.5, 4.8, 7.2, and 20.9 for a p = 2.4, 2.5, 2.7, 
and 2.9, respectively. These values are slightly smaller than in- 
stantaneous jet pow er estimates of M87 being of the order of 
- 10 43 ergs _1 jBicknell & Begelman|[l996l lYoung et alJ 
2002). In general, this demonstrates the ability of future high 
resolution y-ray observations to constrain the energy fraction 
of CRp escaping from the radio plasma. 
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3.4. Radio emissivity of secondary electrons: The 
case of the radio mini-halo in Perseus 

In contrast to y-rays induced by hadronic CRp interactions 
whose spectral shape and normalization is only governed by the 
spectral index a p as free parameter, the resulting radio emission 
from secondary electrons also depends on the morphology and 
strength of the magnetic field B(r). Because only a subsample 
of cooling flow clusters contain radio mini-halos which are not 
outshined by the central AGN we decided to concentrate on 
the Perseus cluster. It has the fortunate property that the radio 
emission due to the central galaxy NGC 1275 is spatially re- 
solved and can be separated from the diffuse emission due to 
the radio-mini halo. 

3.4.1 . Intracluster magnetic fields 

Magnetic fields in galaxy clusters seem to be on the level of 
~ /iG. Indirect estimates of magnetic field strength assuming 
equipartition of energy density of the fields and that of a ra- 
dio synchrotron emitting relativistic electron population give 
low field strengths of ~ 0.1 yuG. Also lower limits on the field 
strength of a comparable level can be derived using the mea- 
surements or upper lim its on IC scattere d CMB phot ons in the 
hard X-ray band JRephaeli et alJfT994t iFusco-Femiano et alJ 
1 19991 lEnfilin et alJ Tl999). Conversely, Faraday rotation mea- 
surements indicate magnetic fields strengths of several fiG in 
typical galaxy c lusters and a few 10/u G in cooling flow re- 
gions of clusters JCarilli & Tavloil2002l for a review). Faraday 
rotation based measurements of the field strength depend on 
estimating the magnetic autocorrelation length from fluctua- 
tions in the Faraday rotation maps. Although the formerly used 
met hods to estimate t his length-scale seem to be questionable 
(EnBlin & Vogt 2003) a refined analysis gives comparable re- 
sults for the magnetic field strengths (Vogt & EnBlin 2003, in 
preparation). 

3.4.2. Comparison of the morphology of radio 
emissivity from secondary electrons 

The radio data was taken from IPedlar etail Jl99(A where we 
neglected the innermost data points because of enhanced con- 
tribution to radio brightness of the radio jet of NGC 1275 and 
the outermost data points due to the limited sensitivity on the 
larger scales of the specific VLA configuration likely leading 
to an artificial decline in the radio surface brightness. The val- 
ues for the azimuthally averaged radio surface brightness were 
converted assuming a two-dimensional Gaussian beam which 
leads to a beam area A bea m = n(4 ln2)^' FWHM A FWHM_ V . 
Figure|4]shows the radial distribution of radio brightness S v (r ± ) 
as a function of impact parameter r± obtained by means of 
eq. J42l> in comparison to the radio data. The CRp adiabatic and 
isobaric model being described in Sect. l3.2l are both shown us- 
ing model parameters of a p = 2.3, Bq = 10 juG, and = 0.5, 
where the latter two parameters refer to eq. (I4U . The normal- 
ization of the radio brightness depends on the assumed scaling 
between CRp and thermal energy density. We fix this scaling 
parameter Xcr p by comparing the simulated radio brightness 



isobaric model 

adiabatic model 




r x [hjl kpc] 

Fig. 4. The radial distribution of radio brightness as a function 
of impact parameter r ± . Shown are the CRp adiabatic and iso- 
baric model for model parameters Bo = 10 yt/G, ag = 0.5, 
and ffp = 2.3 (details are described in the text) as well as 
the azimuthally averaged radio brightness profil e of the the 
Perse us radio-mini halo (data was taken from IPedlar et alJ 
1990). The normalization of the radio brightness depends on 
the assumed scaling between CRp and thermal energy density. 
We fix this scaling parameter Xqr p by comparing the simulated 
radio brightness to the measured data at 24.65 kpc. 

Table 5. Upper limits on the CRp scaling parameter Xcr p in- 
ferred from radio brightness profiles of the radio mini-halo of 
Perseus cluster for different values of Bo, ag, and a p . 



Model 




So foiG] 


a B 


yisobaric 
A CRp 


yadiabatic 
A CRp 


1 


2.3 


10 


0.5 


0.016 


0.006 


2 


2.1 


10 


0.5 


0.014 


0.005 


3 


2.5 


10 


0.5 


0.033 


0.011 


4 


2.7 


10 


0.5 


0.096 


0.031 


5 


2.3 


5 


0.5 


0.027 


0.009 


6 


2.3 


20 


0.5 


0.012 


0.004 


7 


2.3 


10 


0.7 


0.017 


0.006 


8 


2.3 


10 


0.9 


0.019 


0.006 



to the measured data at 24.65 kpc. There is an excellent 
morphological concordance of the isobaric model of CRp and 
the radio data for the radio-mini halo of the Perseus cluster. 
Since the required values of Xqr p are plausible (~ 0.01-0.1, see 
Sect. 13.4.31 . the hadronic secondary CRe model is a very attrac- 
tive explanation for the observed radio mini-halos in cooling 
flow clusters. 

3.4.3. Results on the scaling parameter X C r p using 
radio observations in different models 

By comparing the simulated radio brightness to the measured 
radio data at 24.65 hzl kpc which is the innermost azimuthally 
averaged data point not being outshined by the radio galaxy 
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Fig. 5. The deprqjected scaling parameter Xcr p (/) between 
CRp and thermal energy density in the adiabatic model ap- 
plied to the mini-radio halo of Perseus and presented for mod- 
els which are defined in Table|5] 



cocoon of NGC 1275 we determine the CRp scaling parameter 
Xcr p - Taking this point of reference yields more conservative 
upper limits for Xcr p instead of normalizing by the integrated 
radio surface brightness especially in the case of poorer mor- 
phological matches. The inferred values forXcRp in Tableware 
shown for different combinations of Bo, &b, and a p . 

Deduced values of this scaling parameter Xqr p which are 
obtained by considering only pion decay induced secondary 
electrons resulting from hadronic CRp interactions in the ICM 
reflect upper limits because there are also other mechanisms in 
galaxy clusters leading to relativistic populations of electrons 
(see Sect.^. By analyzing the variations of our model param- 
eters in Table|5]we conclude a weak dependence of Xcr p onaj 
while the magnetic field strength at the cluster center Bo and 
the CRp spectral index a p show a stronger influence on Xcr p - 
The spectral parameter of the magnetic field as impacts mostly 
on the radial extensions of the radio brightness profiles while 
the CRp scaling parameter reflects a degeneracy with respect 
to Z?o and a p . 

Figure[5]shows the scaling parameter XcR P (r) as a function 
of radius r between CRp and thermal energy density in the adia- 
batic model according to eq. J54l > for models defined in Table|5] 
The enhancement of CRp relative to the thermal energy density 
owing to adiabatic compression of the CRp population during 
the formation of the cooling flow can be clearly seen. 

3.5. Constraints derived from the radio halo of Coma 
3.5.1 . Parameter study of the hadronic scenario 

We also applied this formalism of synchrotron radiation emit- 
ted by secondary electrons as presented in Sect. 12.41 to the ra- 
dial distribution of radio brightness in th e radio halo of th e 
Coma cluster using radio data at 1 .4 GHz bv lDeiss et alJ {T997 ). 
Assuming the CRp population to be distributed according to the 
isobaric model, the spatial radio brightness profile obtained by 




100 
r[Kl kpc] 



1000 



Fig. 6. The deprojected CRp scaling parameter XcR P {r) re- 
quired to account for the observed radio halo in Coma within 
the framework of the hadronic scenario. The rising curves with 
increasing radius represent XcR P (r) while the declining curves 
show Xs(r) for the particular choice of a magnetic field be- 
ing fr ozen into the flow and isotropized, i.e. ag — 0.7 (Tribble 
1993). 



this secondary electron model declines too fast with increas- 
ing impact parameter r ± in order to account for the observed 
extended radio halo of Coma. To check whether this shortfall 
of the theoretical model represents a serious problem for the 
hadronic model of radio synchrotron emission we are asking in 
turn for the necessary radial variation of the CRp scaling pa- 
rameter XcR P (r) that is able to explain the observed radio halo. 
Deprojecting the azimuthally averaged observed radio surface 
brightness profile which is described by a /J-model yields (as 
laid down in AppendixIXl 



Mr) 



6/3-1 



2n r c 



(l + P-Irt) 



3/J 



where So = 1.1 mJy arcmkT 2 , r c = • 
By comparing the observed to the theoretically expected ra- 
dio emissivity at each radius we infer the ratio of CRp-to- 
thermal energy density XcR P (r). Figure [6] shows a comparison 
of XcR P (r) and the ratio of magnetic-to-thermal energy den- 
sity Xs(r) = SB(r)/s t h(r) for particular model parameters a p , 
cxb, and Bo- Whereas a p and Bo impact mostly on the nor- 
malization of both scaling parameters XcR P (r) and Xg(r), the 
choice of as governs the relative curvature of these functions: 



'70 



(71) 



kpc, and B = 0.78. 



X B {r) °c « e (r) 



2ar»-l 



is curved in a convex fashion for ag > 0.5 



and exhibits concave curvature for as < 0.5 assuming the clus- 
ter to be isothermal which is a valid approximation for Coma. 
While there are combinations of parameters for which XcR P (r) 
becomes larg er than unity and thus question the hadronic sce- 
nario ( Brunetti 120021) . only small variations in parameter space 
yield plausible values for XcR P (r) (compare Fig.|6j. 

In order to quantify these considerations we perform a pa- 
rameter study to exclude regions of parameter space spanned 
by a p , as, and Bo where the hadronic scenario is challenged 
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to account for the radio halo in Coma. Figure shows the 
resulting contour lines of XcR P (r < 1 hzl Mpc) = 1 and 
^CRp(/ ^ 1 Mpc) = 0.1 in this parameter space. The gra- 
dient of XcR P (r < 1 hjl Mpc) points towards the lower right 
corner in Fig.Qand thus leaves the upper left region of param- 
eter space where the hadronic scenario is able to account for 
the observed radio halo depending on the specific choice of as- 
Since Xg(r < 1 \c\ Mpc) < 0.1 for the entire region of pa- 
rameter space investigated here there are no further constraints 
imposed on the hadronic scenario. 

Choosing the energy density of the magnetic field to de- 
cline like the thermal energy density, i.e. as = 0.5, requires 
XcRp(r) to increase by a factor of less than two orders of mag- 
nitude from the center to the outer parts of the cluster in order to 
reproduce the observed radio halo of Coma. This factor, how- 
ever, is reduced for smaller values of a#. It is further reduced 
due to the non-spherical morphology of Coma, as explained in 
the following. The X-ray emissivity and the radio emissivity 
resulting from hadronic CRp interactions differ in their scaling 
with the electron density according to 

AxO) °c H e O) 2 and (72) 

Mr) <* XcR P (r) « e (r) 2+ff « (I+ V 2) ~ X CRp (r) n e (rf A (73) 

within the framework set by our model and depending on the 
particular choice of as and a p . Thus, any anisotropy like the 
Coma X-ray and radio bridge yields biased azimuthal aver- 
ages when comparing observational to theoretical radio sur- 
face brightness profiles where the latter uses density profiles 
obtained by deprojecting X-ray profiles. Remarkably, this dis- 
crepancy is largest for large values of ag and a p for which we 
infer the tightest limits on the hadronic scenario (cf. Fig. [7Jl 
and thus softens these limits. This results in biased profiles 
of XcR„( r) which increase too strongly towards larger radii 
(cf. also iDolag & E nBlin 2000). Pursuing an approach of av- 
eraging only along the line of sight could attenuate the bias 
JCovoni et all200ll 

An increase of XcR P (r) towards the cluster's periphery is 
indeed observed in cosmological structure formation simula- 
tions due to adiabatic compression inside the cluster which 
increases the thermal pressure at a higher rate than the CRp 
pressure ( Mini ati et all2001a | b|). Bearing in mind th at the CRp- 
to-thermal pressure ratio of lMiniati et alJ J2001albh is obtained 
from volume averages and the energy density stored in mag- 
netic fields declines shallower in comparison to the thermal en- 
ergy density we conclude that our results arising the parameter 
study may be well in agreement with these simulations. 

3.5.2. The spectrum of the Coma radio halo 

One might object that the CRp spectral index should be de- 
termined better owing to radio observations than the range of 
a p - (2, 3) being considered in the previous parameter study 
("Sect. l3~5"TY The following line of argumentation shows, that 
this, on the contrary, is not the case. First, there is an ambigu- 
ity of relating the CRp spectral index a p to the induced syn- 
chrotron spectral index a v which is either a v = a p /2 (Dermer's 
model) or a v = (2a p - l)/3 (fireball model). When compar- 
ing multifrequency observations of diffuse radio emission of 




Fig. 7. Parameter study on the ability of hadronically originat- 
ing CRe to generate the radio halo of Coma. Shown are con- 
tour lines of XcR P (r < 1 tcl Mpc) = 1 (solid) and XcR P (r < 
1 fori Mpc) = 0.1 (dotted) in parameter space spanned by a p 
and Bo for three choices of magnetic field morphology charac- 
terized by as- The contour line of Xcr p = 0.1 for as = 0.3 
has been omitted since it almost coincides with the contour of 
^CRp = 1 for as = 0.7. The lower right corner represents the 
region in parameter space, where the hadronic scenario faces 
challenges for explaining the observed radio halo of Coma. 



the ICM which extends to several GHz the Sunyaev-Zel'dovich 
(SZ) distortion of the spectrum has to be taken care of. At these 
frequencies of the Rayleigh-Jeans part of the Planck spectrum 
the SZ effect amounts to a decrement which introduces a cutoff 
in the r adio spectrum as can be seen in Fig. [8] Following E nfilinl 
(2002) the SZ luminosity reads for Coma in the Rayleigh-Jeans 
part 

1/2 Jy, 



F Coma = _4 j x 10 -3 V 2 



GHz ^70 



(74) 



where vghz = v/GHz. 3 However, the SZ amplitude is uncer- 
tain within a factor of 2 which stems mostly from density pro- 
files being inferred from X-ray observations and extrapolated 
to /?shock- Furthermore, the multifrequency dataset as compiled 
bv lThierbach et all d2003h is inhomogeneous because the solid 
angle over which the observed radio fluxes have been integrated 
may vary among these observations. Finally, the quoted uncer- 
tainties may underestimate the systematic uncertainties which 
e.g. result from incomplete accounting for point source sub- 
traction. 

4. Detectability of y-rays by satellite missions and 
Cerenkov telescopes 

Based on the previous results we discuss the detectability of 
IC emission by secondary CRe and pion decay induced y-ray 
emission by current and future satellite missions as well as op- 
erating and future Cerenkov telescopes. 

3 He r e we corrected for a missing factor of 2 in equation (4) in 
EnBlin 1 2002) and changed the slope of the /3-profile to /? = 0.75 
iBriel et alll992t) . 
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Table 6. Expected limits on the CRp scaling parameter Xcr p by comparing the integrated pion decay induced y-ray flux above 
100 GeV to sensitivity limits of Cerenkov telescopes of T y ,^ p (E > E^) = 10~ 12 y cirT 2 s (E^jlQQ GeV) 1- "'' assuming a 
y-ray spectral index in Dermer's model a y - a p . Note that limits on Xqr p roughly < 0.01 for a p - 2.3 in the isobaric model 
provide good chances to detect y-rays in these particular clusters with new generation Cerenkov telescopes. 



visobaric vadiabatic 

A CRp A CRp 

Cluster o-p = 2.1 a p = 2.3 a p = 2.5 a p = 2.7 a p = 2.1 a p = 2.3 a p = 2.5 a p = 2.7 

A85 3.7 x 1(T 2 6.5 x 1(T 2 2.2x10-' 1.0 2.7 x 10~ 2 4.7 x 10~ 2 1.5x10-' 7.1 x 10~' 

A426 (Perseus) .... 2.5 x 10" 3 4.5 x 10~ 3 1.5 x 10~ 2 7.0 x 10~ 2 2.1 x 10~ 3 3.6 x 10~ 3 1.2 x 10~ 2 5.5 x 10~ 2 

A2199 4.4 x 10~ 2 7.7 x 10~ 2 2.6x10-' 1.2 4.1 x 10~ 2 7.2 x 10~ 2 2.4x10-' 1.1 

A3526 (Centaurus) 1.9 x 10~ 2 3.4 x 10~ 2 1.1 x 10 -1 5.3 x 10~' 1.8x10^ 3.2 x 10~ 2 1.1 x 10 -1 5.0 x 10 _1 

Ophiuchus 4.0 x 10~ 3 7.0 x 10~ 3 2.3 x 10~ 2 1.1 x 10 _1 

Triangulum Australis 1.6 xlO -2 2.8 x 10~ 2 9.3 x 10~ 2 4.4 XlO" 1 

Virgo 4.1 x 10~ 3 7.3 x 10~ 3 2.4 x 10~ 2 1.1 x 10'' 3.8 x 10~ 3 6.7 x 10~ 3 2.2 x 10~ 2 1.0 x 1Q-' 

A1656 (Coma) 7.9 x 10~ 3 1.4xl0~ 2 4.6 x 10~ 2 2.2x10-' 

A2256 4.9 x 10~ 2 8.6xl0~ 2 2.9x10-' 1.4 

A2319 1.5 x 10~ 2 2.6 x 10~ 2 8.9 x 10~ 2 4.2x10-' 

A3571 1.9 xl<T 2 3.4xl0~ 2 1.1 x 10" 1 5.3 X 10" 1 
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Fig. 8. Observed radio halo fluxes of the Coma cluster as com- 

I 1 j 1 

piled by Thierbach et al. (2003). Shown are synchrotron power- 
law spectra for different spectral indices a Y - a p /2. The spec- 
tra are modified at higher frequencies by means of the SZ effect 
while the SZ decrement (a negative flux at the Rayleigh- Jeans 
part of the Planck spectrum) is also shown (dotted). The SZ flux 
is derived from the cluster volume up to the assumed position 
of the accretion shock. 



4.1. Detectability of pion decay induced y-ray and IC 
emission by secondary CRe with INTEGRAL 

The imager IBIS which is the Imager on Board the 
"INTErnational Gamma-Ray Astrophysics Laboratory" 
(INTEGRAL) 4 Satellite covers an energy range from 15 keV 
up to 10 MeV and is capable of high resolution imaging 
(12' FWHM) and source identification. Its spectral sensitivity 
reaches down to 5 x 10~ 8 y s _1 cirT 2 keV -1 (3 cr in 10 6 s, 
AE = E/2) to the continuum at 10 MeV. However, this 



is most probably not sufficient in order to detect the pion 
decay induced y-rays of a particular cluster (compare Fig. [5}. 
Assuming a CRp spectral index of a p = 2.3 and taking the 
results of Table [5] we expect an IC emission of hadronically 
originating CRe in the Perseus cluster of 



^(20keV) = ^F Ic 10- 7 ycm- 
dE 



keV~ 



(75) 



with Tic = 8.4, 4.2, and 2.3 for B = 5 fiG, 10 fiG, and 20 jt/G 
Comparing these results to the post-launch spectral sensitiv- 
ity of 4 x 10~ 6 y s _1 cirT 2 keV -1 to the continuum at 20 keV 
for an observation time of 10 6 s (3 cr detection) there is only 
a minor chance to detect IC emission of CRe. However, for 
steeper spectral indices or a strongly inhomogeneously magne- 
tized environment, IC fluxes can b e enhanced at the expense of 
synchrotron emission according to Enfi lin et al .1 ill 9991) . 



4.2. Possibility of pion decay induced y- 
by GLAST 



y-ray detection 



The "Large Area Telescope" (LAT) onboard the "Gamma- 
ray Large Area Space Telescope" (GLAST) 5 scheduled to 
be launched in 2006 has an angular resolution smaller than 
3.5° at 100 MeV while covering an energy range of 20 MeV 
up to 300 GeV with an energy resolution smaller than 10%. 
Assuming a photon spectral index ofa y =2 for the y-ray back- 
ground the point source sensitivity at high galactic latitude in an 
one year all-sky survey is better than 6x 10~ 9 cirT 2 s _1 for ener- 
gies integrated above 100 MeV. Assuming the radio-mini halo 
in the Perseus cluster mainly to originate from secondary elec- 
trons emitting synchrotron radiation then we expect the CRp 
scaling parameter to be typically one order of magnitude below 
the upper limits obtained by comparing to EGRET data. This 
immediately would imply a good possibility to detect pion de- 
cay induced y-ray emission by GLAST preferentially in nearby 



http://astro.esa.int/Integral/ 



5 http : //glast . gsf c . nasa . gov/science/ 
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cooling flow clusters like Perseus and Virgo. Specifically for 
our secondary model of the radio mini-halo of Perseus, while 
assuming a p - 2.3 in the CRp isobaric model we expect an 
integrated y-ray flux above 100 MeV from Perseus of T y {> 
100MeV)/(ycirT 2 s" 1 ) = 1.3xl0- 8 ,7.4xl0~ 9 , and5.6xl0~ 9 
for Bo = 5yuG, lOjt/G, and 20/iG The expected y-ray flux is 
ever higher when including lower energetic photons. 

4.3. Expected y-ray flux for Cerenkov telescopes 

In the near future there will be different Cerenkov tele- 
scope experiments operating with several telescopes simulta- 
neously and therefore allowing stereoscopic observations. On 
the southern hemisphere there are the "Collaboration between 
Australia and Nippon for a Gamma Ray Observatory in the 
Outback" (CANGAROO) 6 in Australia and the "High Energy 
Stereoscopic System" (HESS) 7 in Namibia. On the north- 
ern hemisphere there will be the "Very Energetic Radiation 
Imaging Telescope Array System" (VERITAS) 8 in Arizona. 
All these telescopes have comparable lower energy thresholds 
of Z?thr = 100 GeV and provide flux sensitivities better than 
T y , ew (E > 100 GeV) = 10" I2 ycrrr 2 s" 1 . On the northern 
hemisphere there will also be the "Major Atmospheric Gamma- 
ray Imaging Cerenkov detector" (MAGIC) 9 on the Canary 
Islands observing with a single dish telescope of 234 m 2 pro- 
viding an even lower energy threshold of E m i n = 30 GeV. 

Following the formalism described in Sect. l3.3.2l and com- 
paring the resulting y-ray flux T y {E > £ t hr) to expected flux 
sensitivities of Cerenkov telescopes T y ^ xp {E > Zs t hrX we ob- 
tain possible upper limits on the CRp scaling parameter Xcr p 
for an integrated volume out to a radial distance of 3 hzl Mpc. 
Table |6] shows constraints for Xcr p using the isobaric and the 
adiabatic model of CRp described in Sect. 13.21 By comparing 
these limits to those obtained by analyzing synchrotron emis- 
sion in the Perseus and Coma cluster (see Table |5J and as- 
suming a substantial contribution of hadronically originating 
CRe to these radio halos there is a realistic chance to detect ex- 
tragalactic pion decay induced y-ray emission in clusters like 
Perseus (A 426), Virgo, Ophiuchus, and Coma (A 1656). 

5. Conclusion 

We investigated hadronic CRp-p interactions in the ICM of 
clusters and simulated the resulting emission mechanisms in 
radio, X-rays, and y-rays assuming spherical symmetry. By 
applying this technique to a sample of prominent clusters of 
galaxies including cooling flow clusters we succeeded in con- 
straining the population of CRp. Especially cooling flow re- 
gions are perfectly suited for constraining non-thermal ICM 
components due to their high gas density and magnetic field 
strength. 

For the first time we developed an analytic formalism to 
describe the 7r°-decay induced y-ray spectrum self-consistently 

6 http : //www . physics . adelaide . edu . au/astrophysics/ 

7 http : //www.mpi-hd.mpg . de/hfm/HESS/HESS .html 

8 http : //veritas . sao . arizona . edu/ 

9 http : //hegral . mppmu . mpg . de/MAGICWeb/ 



for a given differential number density distribution of the CRp 
population being described by a power-law in momentum p p 
and parametrized by the spectral index a p . Assuming a con- 
stant scaling between kinetic CRp energy density and ther- 
mal energy density of the ICM we derived an analytic < T r y~Fy^ 
scaling relation which only applies accurately for isothermal 
clusters. Given the bolometric X-ray luminosity of a particular 
cluster this formula estimates the expected y-ray flux T y ow- 
ing to inelastic cosmic ray ion collisions. From the literature 
we collected electron density and temperature profiles of seven 
cooling flow clusters and four non-cooling flow clusters using 
the Ty-Fx scaling relation to obtain observationally promis- 
ing candidates. We furthermore present formulae describing 
the synchrotron and inverse Compton emission of hadronically 
originating secondary electrons assuming an isotropic distribu- 
tion of magnetic fields following a smooth profile. 

In order to apply this method to our sample of clusters 
of galaxies we introduced three specific models for the spa- 
tial distribution of CRp within cooling flow cluster. In our first 
two scenarios we characterized the kinetic CRp energy density 
fiCRpW to be a constant fraction of the thermal energy den- 
sity ethW of the ICM parametrized by Xcr p . The CRp isobaric 
model assumes the average pressure of CRp not to change dur- 
ing the formation of the cooling flow while the adiabatic model 
hypothesizes this proportionality prior to transition because the 
CRp experience adiabatic compression during the relaxation 
phase. In our third scenario we modeled the resulting distribu- 
tion of CRp diffusion from a central source. By modeling the 
particular y-ray emission of our cluster sample and compar- 
ing to EGRET upper limits we obtained upper bounds on the 
CRp scaling parameter Xcr p = scRp(r)/s±(r). For Perseus and 
Virgo we infer the strongest upper limits which lie in the range 
Xcrp e [0.08,0.18] for different choices of the CRp spectral 
index a p e [2.1, 2.7]. 

Furthermore, the radio emission due to hadronically pro- 
duced secondary electrons emitting synchrotron radiation was 
calculated and resulting radio brightness profiles were com- 
pared to measured data of the radio-mini halo of Perseus as 
well as the radio halo of Coma. In the case of Coma our CRp 
profiles characterized by a flat CRp scaling parameter Xcr p are 
not able to reproduce the observed radio profiles particularly 
in the peripheral regions of the cluster. In the following we ad- 
justed the radial behavior of XcR P (r) such that the synchrotron 
emission resulting from hadronic CRe is able to account for 
the observed radio surface brightness profile and thus allow- 
ing for an additional degree of freedom. The resulting increase 
of XcR P (r) for larger radii could be due to adiabatic compres- 
sion which increases the thermal energy density at a higher 
rate than the CRp energy density. Even more important, the 
aspherical Coma cluster morphology reduces the required ra- 
dial increase in XcR P (r). By exploring the accessible parameter 
space spanned by parameters describing the magnetic field and 
the spectral index of the CRp population we identify regions 
where the hadronic scenario is able to reproduce the observed 
radio profiles preferentially for an energy density of the mag- 
netic field which declines shallower than the thermal energy 
density. We conclude that the secondary model for radio halos 
is still viable. 
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In the case of the Perseus mini-radio halo, we conclude up- 
per limits on Xqr p which are ranging for the isobaric model 
of CRp within the interval Xqr ? e [0.01,0.1] for conserva- 
tive combinations of values of the magnetic field B and the 
CRp spectral index a p while upper limits for the CRp adiabatic 
model are typically half an order of magnitude below. By com- 
paring calculated radio brightness profiles to measured data of 
the radio-mini halo in Perseus, we found excellent morpholog- 
ical agreement between the CRp isobaric model and the radio 
data especially for the choice of Bq - 10 fiG, as = 0.5, and 
a v — 2.3. In the course of this paper we argued that this specific 
choice of parameters for the magnetic fields in cooling flow 
clusters is also preferred by experiments like Faraday rotation 
measurements and cosmological cluster simulations including 
magnetic fields. A discussion of different acceleration mecha- 
nisms of CRp such as structure formation shocks, supernovae 
remnants, and injection by active radio galaxies supports also 
a value of a v close to the inferred one. Because of the required 
moderate CRp energy density we propose synchrotron radia- 
tion by non-thermal secondary electrons from hadronic inter- 
actions as a likely explanation of radio mini-halos. In order to 
scrutinize this model we provide predictions of y-ray fluxes for 
Cerenkov telescopes as well as the INTEGRAL and GLAST 
satellites. 

Finally, we analyzed the possibility of detecting such pion 
decay induced y-ray and IC emission by current and future 
satellite missions as well as new generation Cerenkov tele- 
scopes. Depending on the CRp spectral index, the fragmenta- 
tion of the spatial distribution of the magnetic field as well as 
its field strength, it will be difficult for INTEGRAL to detect 
the IC emission of the hadronically originating secondary CRe 
while GLAST has the potentiality to detect the distinct signa- 
ture of the pion decay induced y-ray emission preferentially 
in nearby cooling flow clusters. By investigating the opportu- 
nity of detecting extragalactic y-rays by Cerenkov telescopes 
we argued in favor of four candidate clusters (Perseus (A 426), 
Virgo, Ophiuchus, and Coma (A 1656)) which are especially 
suited to detect hadronically originating y-ray emission. 
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Appendix A: Deprojection of X-ray surface 
brightness profiles represented by double 
/3-models 

Owing to the enhanced electron density in the central region the 
X-ray surface brightness profile Sx(r±) in cooling flow cluster 
can be represented by double B models, 

Sf+l/2 



Sx(rJ = £ S* 



1 + 1^ 



(A.l) 



where the X-ray surface brightness profile is a line of sight pro- 
jection of the squared electron density and the cooling function 
relative to the squared electron density Ax(T e ), 



Sx(rJ 



+ Z A; 



rn 2 e (r)A x [T e (r)] 



(A.2) 
(A.3) 



r 2 - r 2 



Thus the electron density n e (r) can be derived from Sx(r±) by 
inverting the Abel equation 



n 2 Jr)A x [T e (r)] 



i d r 

n r dr J r 

i r 

n Jr 



dy 



ySx(y) 

yjy2 _ r 2 



y/yl - f-2 



(A.4) 
(A.5) 



where the prime denotes the derivative. For the second equation 
we used that n e (r) is bounded for r — > oo. Using eq. ( IA.U this 
equation can be solved analytically yielding 



K(r) = 7 



1 



s 



Si 



6B t - 1 



A x LT e (r)]^2,rr C( ( 1+r 2 /r 2) 



i,3fl-|, (A.6) 



where S(a, b) deno tes the beta-function 

dAbramowitz & Steeunl 1 1965b . Provided the central den- 
sity « e (0) is known and assuming furthermore the special case 
of equality of the two B parameter, B\ — B2, we arrive at the 
following compact formula for the electron density profile 
n e (r) 



n e (r) 



A x rr e (0)] 



AxLT e (r)] 
Hi = n e (0) 



( 2 c 



(=i 

s-l/2 



2\-3>S 



1/2 



S;r c 



(A.7) 



(A.8) 



Generalizing to n-fold y6-profiles can be obtained by means of 
induction. 
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